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Abstract 

Isoscalar and isovector monopole oscillations that correspond to giant resonances 
in spherical nuclei are described in the framework of time-dependent relativistic 
mean-field theory. Time-dependent and self-consistent calculations that reproduce 
experimental data on monopole resonances in ^ospi^ show that the motion of the 
collective coordinate is regular for isoscalar oscillations, and that it becomes chaotic 
when initial conditions correspond to the isovector mode. Regular collective dy- 
namics coexists with chaotic oscillations on the microscopic level. Time histories, 
Fourier spectra, state-space plots, Poincare sections, autocorrelation functions, and 
Lyapunov exponents are used to characterize the nonlinear system and to identify 
chaotic oscillations. Analogous considerations apply to higher multipolarities. 
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1 Introduction 



In the last two decades many studies have been reported in which the atomic nucleus 
has been used as a laboratory, both experimentally and theoretically, for investigating 
the transition from order to chaos in quantum dynamical systems (for a recent review 
see [|l|). Most of these studies have concentrated on two major aspects: (i) generic 
signatures of chaos in local fluctuations and correlations of nuclear level distributions, and 
(ii) chaos in microscopic and collective dynamics of realistic many-body systems. In the 
first case, signatures of quantum chaos have also been studied in the complicated structure 
of wave functions and randomness of matrix elements of physical operators. On the other 
hand, the nature of collective nuclear dynamics has been investigated with particular 
emphasis on the stability of low-lying nuclear modes in relation to one-body dissipation 
caused by the deformation of the nuclear potential, cluster effects and Coriolis forces. 
The most elementary collective modes in nuclei are giant resonances. These are highly 
collective nuclear excitations in which a large fraction of nucleons participate. They can 
be described as damped harmonic/anharmonic density oscillations around the equilibrium 
ground-state. Giant resonances occur over the whole periodic table and their characteristic 
parameters are smooth functions of the mass number. A mean-field model therefore 
provides an appropriate framework for the description of giant resonances. Regular and 
chaotic dynamics in giant nuclear oscillations has been the subject of a number of studies. 
What has emerged as a very interesting result is that an undamped collective mode may 
coexist with chaotic single-particle motion. It appears that the slowly vibrating self- 
consistent mean field created by the nucleons averages out the random components in 
their motion. In all investigations the motion of only one type of particles has been 
considered. That is, only the dynamics of isoscalar collective modes. 
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In the present article we study the difference in the dynamics of isoscalar and isovec- 
tor collective modes. In particular, we describe isoscalar and isovector monopole oscil- 
lations in spherical nuclei, but analogous considerations apply to higher multipolarities. 
The dynamics of nuclear vibrations is analyzed in the framework of time-dependent rela- 
tivistic mean- field theory (TD RMFT). The model represents a relativistic generalization 
of the time-dependent Hartree-Fock approach. Nuclear dynamics is described by the si- 
multaneous evolution of A single particle Dirac spinors in the time- dependent mean fields. 
Frequencies of eigenmodes are determined from a Fourier analysis of dynamical quantities. 
In this microscopic description, self-consistent mean-field calculations are performed for 
static ground-state properties, and time-dependent calculations for monopole excitations. 
Because the time-evolution is calculated self-consistently, the system is intrinsically non- 
linear. An advantage of the time-dependent approach is that no assumption about the 
nature of the mode of vibrations has to be made. 

The article is organized as follows. In Sec. 2 we present the essential features of the 
time-dependent relativistic mean-field model, as well as some details of its application 
to spherical nuclei. Time-dependent calculations of isoscalar and isovector monopole 
oscillations in ^^O and ^''^Pb are described in Sec. 3. Results of a number of diagnostic 
tests that are used to identify chaotic oscillations in the nuclear system are discussed and 
compared with earlier work. A summary of our results is presented in Sec. 4. 

2 Outline of the model 

The dynamics of collective vibrations in spherical nuclei will be described in the framework 
of relativistic mean-field theory [@, @, Hj. In relativistic quantum hadrodynamics the 
nucleus is described as a system of Dirac nucleons which interact through the exchange 
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of virtual mesons and photons. The Lagrangian density of the model is 

The Dirac spinor denotes the nucleon with mass m. m^, m^, and rrip are the masses 
of the cr-meson, the w-meson, and the p- meson, and g^, Qui, and Qp are the corresponding 
coupling constants for the mesons to the nucleon. U{a) denotes the nonlinear a self- 
interaction 

U{a) = \mla' + ^g,a' + \g,a\ (2) 

and fi''^, R''^, and F'^^ are field tensors 0. 

The coupled equations of motion are derived from the Lagrangian density (||). The 
Dirac equation for the nucleons: 



idti/ji 



a I -iV - g^LJ - gpTp - e^^— j + P{m + g^a) 



+gu^^o + QpTpQ + e Aq 



(3) 



and the Klein-Gordon equations for the mesons: 



(a2-A + m2)a = - (72^' - (73^' (4) 

{d^^ - A + ml)u^, = g^j^ (5) 

{d^-A + ml)p, = gpj, (6) 

(5'-A)^m = ej^. (7) 

In the relativistic mean field approximation the A nucleons described by single-particle 
spinors ipi, = 1) 2, A), are assumed to form the A-particle Slater determinant |$), 



and to move independently in the classical meson fields. The sources of the fields, i.e. 
densities and currents, are calculated in the no-sea approximation : - the scalar density 

A 

Ps = ^AAi (8) 

i=l 



the isoscalar baryon current 



the isovector baryon current 



3' = jl^^r^^1 (9) 
j=l 

= j2i'il^r^^1 (10) 

i=l 

- the electromagnetic current for the photon-field 

3^ra = E^.7^^^^„ (11) 
i=l ^ 

where the summation is over all occupied states in the Slater determinant |$). Negative- 
energy states do not contribute to the densities in the no-sea approximation of the sta- 
tionary solutions. However, negative energy contributions are included implicitly in the 
time-dependent calculation, since the Dirac equation is solved at each step in time for a 
different basis set [§], 0] . Negative energy components with respect to the original ground 
state basis are taken into account automatically, even if at each time the no-sea ap- 
proximation is applied. It is also assumed that nucleon single-particle states do not mix 
isospin. 

The ground state of a nucleus is described by the stationary self-consistent solution 
of the coupled system of equations (0)-(0), for a given number of nucleons and a set of 
coupling constants and masses. The solution for the ground state specifies part of the 
initial conditions for the time- dependent problem. For a given set of initial conditions, 
i.e. initial values for the densities and currents in Eqs. (|8|-[ll|), the model describes the 



time evolution of A single-particle wave functions in the time-dependent mean fields. 



Retardation effects for tlie meson fields are not included, i.e. the time derivatives d'f 
in the equations of motions for the meson fields are neglected. This is justified by the 
large masses in the meson propagators, causing a short range of the corresponding meson 
exchange forces. Since there is no systematic procedure to derive the initial conditions 
that characterize the motion of a specific mode of the nuclear system, the description of 
nuclear dynamics as a time-dependent initial-value problem is intrinsically semi-classical. 
The theory can be quantized by the requirement that there exist time-periodic solutions of 
the equations of motion, which give integer multiples of Planck's constant for the classical 
action along one period |Q. For giant resonances the time-dependence of collective 
dynamical quantities is actually not periodic, since generally giant resonances are not 
stationary states of the mean-field Hamiltonian. The coupling of the mean-field to the 
particle continuum allows for the decay of giant resonances by direct escape of particles. In 
the limit of small amplitude oscillations, however, the energy obtained from the frequency 
of the oscillation coincides with the excitation energy of the collective state. In Refs. 0, 



TT[| we have shown that the model reproduces reasonably well the experimental data 



on giant resonances in spherical nuclei. 

In the present article we apply the model to isoscalar and isovector monopole os- 
cillations in spherical nuclei. In this microscopic description, self-consistent mean-field 
calculations are performed for static ground-state properties, and time-dependent cal- 
culations for monopole excitations. Because of the self-consistent time-evolution, the 
system is intrinsically non-linear. For a system with spherical symmetry, the nucleon 
single-particle spinor is characterized by the angular momentum j, its z-projection m, 
the parity vr, and the isospin ts = ±| for neutrons and protons 

^(r,M,t3) = - { -f^.^'^H'^'l ] e-^^^xAh). (12) 
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Xt is the isospin function, the orbital angular momenta I and / are determined by j and 
the parity vr, /(r) and g{r) are radial functions, and ^ijm is the tensor product of the 
orbital and spin functions 

^lj^{e,ip,s) = <^ ^ ^ > yimii^^V) Xmsis). (13) 

nisTni ^ 

In order to excite monopole oscillations, the spherical solution for the ground-state has 
to be initially compressed or radially expanded by scaling the radial coordinate. The 
amplitudes and g^"^ of the Dirac spinor are defined 

n^'ir^on) = (1 + a) /(r), g^°^{r^,.) = (1 + a) g{r) . (14) 

The new coordinates are 

rmon = (1 + a) r. (15) 

For isoscalar oscillations the monopole deformations of the proton and neutron densities 
have the same sign. To excite isovector oscillations, the initial monopole deformation pa- 
rameters of protons and neutrons must have opposite signs. After the initial deformation 



(|l^), the proton and neutron densities have to be normalized. 

The time- dependent Dirac equation (^ is reduced to a set of coupled first-order 
partial differential equations for the complex amplitudes / and g of proton and neutron 
states 

^^tf = {Vo + g.a)f + idr---tVr)g (16) 

r 

idtg = {Vo - ga(r -2m)g - {dr + - -iVr)f, (17) 

r 

where k, = ±(j + |) for j = / =F |, and the indices and r denote the time and radial 
components of the vector field 

= g^LU^^ + g.pt^T, + e^i^A^ (18) 



For a given set of initial conditions, the equations of motion propagate the nuclear system 
in time. The potentials are solutions of the Klein-Gordon equations 

+ = ^dr), (19) 

are meson masses for (p = a, u and p, and zero for the photon. The source terms 



are calculated from (|)-([TT|) using in each time-step the latest values for the nucleon 
amplitudes. The meson fields and electromagnetic potentials are calculated from 



oo 

G4r,r')s4ry^dr' , (20) 
where for massive fields 

C (r,r') = -^i^(e-^l'--'l_e— (21) 
2m^ rr ^ ' 

and for the Coulomb field 

Gc{r, r') = - for r > r' 
r 

Gc{r,r') = ^ for r < r'. (22) 

The collective dynamical variables that characterize vibrations of a nucleus are defined as 
expectation values of single-particle operators in the time-dependent Slater determinant 
|<l>(t)) of occupied states. In the framework of the TDRMF model the wave-function of 
the nuclear system is a Slater determinant at all times. For isoscalar monopole vibrations, 
the time-dependent monopole moment is defined: 

{r\t)) = jimir'Mt)). (23) 

The corresponding isovector monopole moment is simply 

< ^p(^) > - < rl{t) > . 

Fourier transforms of the collective dynamical variable determine the frequencies of eigen- 
modes. 
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3 Isoscalar and isovector monopole oscillations 

The study of isoscalar monopole resonances in nuclei provides an important source of 
information on the nuclear matter compression modulus -ft'nm- This quantity is crucial in 
the description of properties of nuclei, supernovae explosions, neutron stars, and heavy 
ion collisions. Modern non-relativistic Hartree-Fock plus RPA calculations, using both 
Skyrme and Gogny effective interactions, indicate that the value of should be in the 



range 210-220 MeV |12]. In the framework of relativistic mean field theory on the other 
hand, calculations based on an effective interaction [|T^ with nuclear matter compression 
modulus -f^nm ~ 250 — 270 MeV are in better agreement with available data on spherical 
nuclei |]TT|. The complete experimental data set on isoscalar giant monopole resonances 
(GMR) has been recently analyzed by Shlomo and Youngblood [0. In Fig. 1 we display 
results of time-dependent relativistic mean-field calculations of isoscalar and isovector 
oscillations in ^°^Pb. The experimental isoscalar GMR energy in ^osp]^ jg ^gjj established 
at 13.7 ± 0.3 MeV. Experimental data on isovector giant monopole resonances are much 
less known. The systematics of excitation energies does not, in general, depend on the 
nuclear matter compression modulus. The experimental value for the isovector GMR 
in 208pb: 26 ± 3 MeV [|14[. We have calculated the ground-state with the NLl p 



parameter set of the effective Lagrangian. This effective force has been extensively used 
in the description of properties of finite nuclei along the valley of beta stability P . For 
NLl (Knin = 211.7 MeV), we expect the calculated excitation energy for the isoscalar 
mode to be approximately 1-2 MeV lower than the average experimental value. However, 
the precise value is not important in the present consideration. For the initial deformation 
parameter in (|T^ we have used a = 0.2. In the isoscalar case both proton and neutron 
densities are radially expanded, while for the isovector mode the proton density is initially 
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compressed by the same amount. Therefore although the initial conditions are different, 
in both cases we follow the time evolution of the same system. 

In Fig. la we plot the time history of the isoscalar monopole moment (r^(t)), and 
in Fig. lb the corresponding isovector moment is shown. The isoscalar mode displays 
regular undamped oscillations, while for the isovector mode we observe strongly damped 
anharmonic oscillations. On the right hand panels we plot the corresponding Fourier 
power spectra. The evolution time Tfinai = 3000 fm/c determines the numerical resolution 
AE = 2 tt/ic/ Tfinai =~ 0.4 MeV in the frequency domain. The Fourier power spectra 
are calculated by FFT using data windowing with the Welch window function. We have 
verified that the positions of the main peaks do not depend on the choice of the window 
function. We estimate the numerical accuracy to ~ 0.5 MeV. As one would expect for a 
heavy nucleus, there is very little spectral fragmentation in the isoscalar channel, and a 
single mode dominates at the excitation energy of 11 MeV. The Fourier spectrum of the 
isovector mode is strongly fragmented. However, the main peaks are found in the energy 
region 25 — 30 MeV, in agreement with the experimental data. 

For the isoscalar mode, the time history of the monopole moment and the Fourier 
spectrum show that the oscillations of the collective coordinate are regular. On the 
other hand, the appearance of a broad spectrum of frequencies seems to indicate that the 
isovector oscillations are chaotic. A diagnosis of chaotic vibrations would imply that one 
has a clear definition of such motion. For a quantum system however, the concept of chaos, 
especially in time-dependent problems, is not well defined. And although our description 
of nuclear vibrations is semi-classical, quantum effects like the Pauli principle are present 
in the initial conditions and during the dynamical evolution. A number of diagnostic tests 
can help to identify chaotic oscillations in physical systems [jl6|, and some of them can be 
applied in the present consideration. In Figs. 2 - 4 we display some additional qualitative 



10 



measures which can be used to characterize the response of our nonhnear system. In Fig. 
2 we have constructed the two-dimensional time-delayed pseudo-phase space for isoscalar 
(a), and isovector (b) oscillations shown in Fig. 1. Since information is available on the 
time evolution of just one variable, the collective coordinate, one plots the signal versus 
itself, but delayed or advanced by a fixed time constant [< r^(i) >, < r'^{t+T) >]. The idea 
is that the signal < r'^{t + T) > is related to < r^(i) >, and should have properties similar 
to those in the classical phase plane [< r^(t) >, < P{t) >]. The choice of the time delay r 
is not crucial, except to avoid a natural period of the system. If the motion is chaotic, the 
trajectories in the phase space do not close. For the pseudo-phase planes shown in Fig. 2, 
we have taken r = 20 fm/c. The phase space trajectories for the isoscalar mode are closed 
ellipses, indicating regular oscillations. For the isovector oscillations on the other hand, the 
trajectories are completely chaotic. The strong damping results from one-body processes: 
(i) escape of nucleons into the continuum states and (ii) collisions of the nucleons with the 
moving wall of the nuclear potential generated by the self-consistent mean fields. In Fig. 3 
we display the corresponding Poincare sections constructed from three-dimensional time- 
delayed pseudo-phase space [< r^(t) >, < r^(t + r) >, < r^(t + 2r) >], with r = 20 fm/c. 
The Poincare maps are shown in the planes: < r^(t) >= 35 frri^ for the isoscalar mode, 
and < r'^{t) >— —9.55 fm^ for the isovector mode, respectively. The Poincare map for 
the isoscalar mode consists of two sets of closely located points and therefore confirms 
regular oscillations. For the isovector oscillations the Poincare map appears as a cloud of 
unorganized points in the phase plane. Such a map indicates stochastic motion. Another 
measure that is related to the Fourier transform is the autocorrelation function 



When the signal is chaotic, information about its past origins is lost. This means that 




(24) 
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A{t) ^ as r — * oo, or the signal is only correlated with its recent past. It is also expected 
A{t) of a chaotically modulated signal to be an irregularly modulated waveform. The 
autocorrelation functions for isoscalar and isovector oscillations are shown in Fig. 4. The 
normalization is A[t = 0) = 1. For the isovector mode A{t) shows a rapid decrease and 
the envelope appears as an irregular waveform. Therefore also this quantity indicates that 
the dynamical variable displays chaotic oscillations for the isovector mode. 

In our first example we have considered monopole oscillations in ^''^Pb, for which 
there exist experimental data on both isoscalar and isovector giant monopole resonances. 
Thus our description of the dynamics of monopole oscillations is not just a completely 
artificial model, but in fact corresponds to an experimentally observed physical system. 
On the other hand ^osp]^ jg large and complicated system, in which many single-nucleon 
orbitals contribute to the collective coordinate. In order to study in more detail the dy- 
namics, it is more convenient to consider a light nucleus: ^^O. In Fig. 5 we compare 
the isoscalar and isovector monopole moments and the associated Fourier power spectra. 
For the isoscalar mode a modulation of the signal is observed, and two strong peaks at 
approximately 20 MeV excitation energy are found in the Fourier spectrum. A fragmen- 
tation in the Fourier spectrum is expected for such a light nucleus. The behavior of the 
collective coordinate for the isovector mode appears to be chaotic. The Fourier spectrum 
is strongly fragmented in the energy region 15 — 45 MeV. The pseudo-phase spaces are 
compared in Fig. 6. The results are similar to those for ^osp]^ g^j^^ indicate that isovector 
oscillations are chaotic. 

In Refs. |T7| , p!8|] Blocki et al. analyzed the behavior of a purely classical gas of 
noninteracting point particles enclosed in a hard-wall container, which undergoes peri- 
odic and adiabatic shape oscillations. The wall motion is not modified by the collisions 
with the particles, i.e. the dynamics is not self- consistent. Analyses of Poincare sections 
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and of Lyapunov exponents showed that for low deformations (/ = 2) the motion of par- 
ticles is regular, whereas for higher multipolarities (/ > 3) the scattering of segments 
of wall with positive curvature leads to divergence between trajectories and therefore to 
chaotic motion. A somewhat different result, especially if applied to the nuclear system. 



was found in Ref. |19|. Conditions under which nucleons inside a nucleus can undergo 
chaotic motion were studied. A self-consistent model of separable forces was used in a 
semi classical approximation of the time-dependent Hartree-Fock equation. The test par- 
ticle method was used to solve the Vlasov equation for the time-evolution of the density 
matrix. Isoscalar quadrupole and octupole oscillations were investigated. It was shown 
that, both for quadrupole and octupole deformations, chaotic single-particle dynamics 
leads to regular motion of the collective coordinate, the multipole moment of deforma- 
tion. The origin of the chaotic single-particle dynamics was attributed to the exchange 
of energy between the motion of the individual test particles and the collective motion of 
the multipole coordinate. In particular, it was stressed that self-consistency is essential 
for the generation of regular dynamics from an ensemble of single-nucleons with chaotic 
trajectories. The relationship between chaoticity at the microscopic level and dissipation 



of the collective degrees of freedom was also investigated in Ref. |2y]. Classical parti- 
cles were confined in a two-dimensional nuclear potential whose walls undergo periodic 
shape oscillations of fixed multipolarity. The collective variable appears explicitly in the 
Hamiltonian as an additional degree of freedom. A fully self-consistent description of 
the dynamics of particle motion and the collective coordinate was employed. Isoscalar 
monopole oscillations were studied and it was shown that self-consistency induces chaotic 
single-particle motion. In particular, it was demonstrated how the coupling between the 
collective variable and single-particle dynamics induces macroscopic dissipation. 

In a quantum description of nuclear dynamics, such as the one based on time- 
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dependent rclativistic mean-field theory, we cannot follow the trajectories of individual 
nucleons. In order to better understand the difference between the isoscalar and isovector 
modes at the microscopic level, we consider the individual contributions of single-nucleon 
orbitals to the collective coordinate, the time-dependent monopole moment. In the case of 
1^0 {Z — N — 8), protons and neutrons occupy only three spherical levels in the ground 
state: Is 1/2, Ip 3/2, and Ip 1/2. Monopole oscillations do not break the degeneracy, i.e. 
spherical symmetry is conserved during the time evolution of the nuclear system and we 
can follow the dynamics of each level. In Fig. 7 we display the time histories of the expec- 
tation values of for the three neutron levels. In Fig. 7a we plot < r'^{t) >si/2,p3/2,pi/2 
for the isoscalar mode, for which the oscillations of the collective coordinate are shown in 
Fig. 5a. Fig. 7b corresponds to the isovector monopole oscillations displayed in Fig. 5b. 

In contrast to the collective isoscalar and isovector moments shown in Fig. 5, in Fig. 
7a and 7b we plot the same quantities, expectation values of for the neutron single- 
particle orbitals. Also the initial conditions for neutron motion are the same in both cases. 
The only difference is the initial condition for protons. In the isoscalar case protons os- 
cillate in phase with neutrons. Proton and neutron oscillations are out of phase for the 
isovector mode. In Fig. 7a one observes regular modulated oscillations for all three neu- 
tron levels. For the isovector case the time histories display chaotic oscillations. In Fig. 8 
we display the corresponding Fourier power spectra and the largest Lyapunov exponents. 
Comparing the Fourier spectra, we notice that the differences between isoscalar and isovec- 
tor oscillations are especially pronounced for the Is 1/2 level. Much more fragmentation 
is observed in the isovector spectrum, and the strength is shifted toward higher frequen- 
cies. The Fourier spectra for isovector oscillations display pronounced peaks in the higher 
frequency region. From the time series of the expectation values < r^{t) >si/2,p3/2,pi/2, 
we have also calculated the largest Lyapunov exponents, shown in the right hand panels 
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of Fig. 8. Lyapunov exponents provide a qualitative and quantitative characterization of 
dynamical behavior. They are related to the exponentially fast divergence or convergence 
of nearby orbits in phase space. A system with one or more positive Lyapunov expo- 
nents is defined to be chaotic. The magnitude of the exponents reflects the time scale 
on which system dynamics becomes unpredictable. The largest exponents are displayed 
as function of evolution time in Fig. 9. They have been calculated by the method of 



Ref. ||2T[, which allows the estimation of non-negative Lyapunov exponents from a time 
series. For the Is 1/2 level the calculated exponent for oscillations that correspond to the 
isoscalar collective mode is consistent with the value zero, indicating regular motion. For 
the isovector mode, the Lyapunov exponent is positive and large. Therefore, oscillations 
of < r^(t) >si/2 are chaotic for collective isovector vibrations. For the Ip 3/2 and Ip 1/2 
levels all calculated exponents are positive. For the isoscalar mode the small positive 
values reflect the observed fragmentation in the Fourier spectrum. Those that correspond 
to isovector collective mode are much larger and indicate that the dynamics is indeed 
chaotic. The corresponding Fourier spectra and largest Lyapunov exponents for proton 
levels are very similar. The only difference is that, in addition to meson exchange forces, 
the protons also feel the long range Coulomb interaction. 

Similar results are obtained for ^°^Pb. We have calculated the Fourier spectra and 
largest Lyapunov exponents for the time series of expectation values < > for the 22 
single-particle spherical shell-model neutron states. Typical results are shown in Fig. 
9 for the orbitals Ih 11/2, Ih 9/2, and 2/ 7/2. The column on the left corresponds to 
isoscalar monopole collective oscillations, on the central panels the results for the isovector 
collective mode are shown. Both for isoscalar and isovector oscillations (Fig. 1), the 
initial conditions for neutron motion are the same. All calculated Lyapunov exponents 
are positive, and therefore seem to indicate that the underlying microscopic dynamics 
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is chaotic. In general, we observe much more fragmentation in the Fourier spectra of 
oscillations that correspond to the isovector mode, and the resulting Lyapunov exponents 
have also larger values. The interesting result of course is that for the isoscalar case we 
observe regular oscillations of the collective variable. This would be in agreement with the 
results reported in Ref. [|19| , where chaotic single-particle motion was found in coexistence 



with regular collective dynamics. For the isovector mode on the other hand, the collective 
monopole moment displays chaotic oscillations. This is explained by the fact that protons 
and neutrons effectively move in two self-consistent potentials that oscillate out of phase. 
For example, when neutrons move inward, they scatter on the potential wall with positive 
curvature that is created by protons moving outward. This will lead to pseudo-random 
motion of the nucleons and dissipation of collective oscillations. 

4 Conclusions 

Isoscalar and isovector collective monopole oscillations in spherical nuclei have been an- 
alyzed within the framework of time-dependent relativistic mean field theory. These 
oscillations correspond to the most elementary collective nuclear modes, the giant res- 
onances. In order to investigate the dynamics of collective vibrations, we have ana- 
lyzed time-dependent and self-consistent calculations that reproduce the experimental 
data on monopole giant resonances in spherical nuclei. Because of the self-consistent 
time-evolution, the nuclear system is non-linear and one expects chaotic dynamics for 
specific initial conditions. In particular, we have studied the difference in the dynamics 
of isoscalar and isovector collective modes. Time histories, Fourier spectra, state-space 
plots, Poincare sections, autocorrelation functions, and Lyapunov exponents have been 
used to characterize the nonlinear system and to identify chaotic oscillations. It has been 
shown that the oscillations of the collective coordinate can be characterized as regular 
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for the isoscalar mode, and that they become chaotic when initial conditions correspond 
to the isovector mode. Our results also confirm the conclusions of a number of studies, 
which have shown how a regular collective mode can coexist with chaotic single-particle 
dynamics. However, we have shown that this is the case only for isoscalar modes, that 
is, only if one considers the motion of a single type of particles. When protons and neu- 
trons move out of phase, as it happens for isovector modes, the resulting dynamics of the 
collective coordinate exhibit chaotic behavior. Of course, analogous considerations apply 
also for higher multipolarities, for example isoscalar and isovector giant quadrupole reso- 
nances. Due to numerical problems, for quadrupole oscillations our computer codes could 
not propagate the nuclear system to very long times. However, a preliminary analysis of 
isoscalar and isovector quadrupole oscillations in ^^O ( Tfi^ai — 500 fm/c) have shown 
that results similar to the monopole resonances could be expected. Work along these lines 
is in progress. 
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Figure Captions 

Fig. 1 Time- dependent isoscalar < > (a) and isovector < Tp > — < > (b) monopole 
moments and the corresponding Fourier power spectra for ^°®Pb. The parameter 
set of the effective Lagrangian is NLl. 

Fig. 2 Two-dimensional pseudo-phase space for isoscalar (a) and isovector (b) monopole 
oscillations in ^°^Pb. 

Fig. 3 Poincare sections for isoscalar (a) and isovector (b) monopole oscillations in ^°^Pb. 

Fig. 4 Autocorrelation functions for isoscalar (a) and isovector (b) monopole oscillations 
in 208pb. 

Fig. 5 Time-dependent isoscalar (a) and isovector (b) monopole moments and the Fourier 
spectra for ^^O. The effective force is NLl. 

Fig. 6 Two-dimensional pseudo-phase space for isoscalar (a) and isovector (b) monopole 
oscillations in ^^O. 

Fig. 7 Time-dependent expectation values < > for the three single neutron spherical 
shell model states Is 1/2, Ip 3/2, and Ip 1/2 in ^^O. (a) corresponds to the isoscalar 
oscillations of the collective coordinate. In panel (b) results for the isovector collec- 
tive mode are shown. 

Fig. 8 Fourier spectra and the largest Lyapunov exponents for the time-dependent ex- 
pectation values < > in the three single neutron spherical shell model states 
Is 1/2, Ip 3/2, and Ip 1/2 in ^^O. Dashed hues indicate the isoscalar collective 
monopole mode, solid lines correspond to isovector collective oscillations. 
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Fig. 9 Fourier spectra and the largest Lyapunov exponents for the time-dependent expec- 
tation values < > in the Ih 11/2, Ih 9/2, and 2/ 7/2 single neutron spherical shell 
model states in ^°^Pb. Results that correspond to the isoscalar collective monopole 
mode are shown on the left, those for the isovector mode are in the center. Lyapunov 
exponents are displayed on right hand side panels. 
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